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Abstract 

The theoretical description of the heterogeneous nucleation kinet- 
ics is presented. This description takes into account the perturbation 
of the vapor phase initiated by the growing droplets. The form of the 
density profile around the growing droplet is analyzed which leads to 
some special approximations. Then the process of nucleation in the 
whole system is described. As the result all main characteristics of 
the process are determined analytically. 

Among numerous examples of the first order phase transitions of 
condensation is extracted by a relative simplicity. This case is well investi- 
gated experimentally and is traditionally regarded as the base for application 
of some new theoretical methods. A "classical theory of condensation" (see, 
for example 0) gives the solid ground for further theoretical constructions. 
Numerous modifications and reconsiderations (see, for example H) allow to 
consider an example of condensation as well analyzed both theoretically and 
experimentally. 

One has to stress that practically all investigations were intended to de- 
termine the rate of nucleation and didn't give the global picture of the phase 
transition. Some theoretical descriptions of a global evolution appeared later 
than the classical theory of nucleation has been created and they were not 
so numerous. Among them one can extract the descriptions of a metastable 
phase decay by Wakeshima [Q, Segal' and Kuni, Grinin [[?[]. A process of 
condensation which occurs during a smooth variation of external conditions 
was considered in ||. Nevertheless the theoretical description of the global 
picture of a whole condensation kinetics have ignored some important fea- 
tures of this process. Namely, an exhaustion of a metastable phase near a 
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growing embryo of a new phase have not been taken into account in a proper 
way. Certainly, this exhaustion is explicitly taken into account in expression 
for the rate of embryos growth in continuous model (i.e. in the "diffusion 
regime" of the embryos growth). This effect was carefully analyzed both in 
stationary and non- stationary aspects in many publications mainly in the 
field of mechanics of continuous media. But the presence of the gap of a 
metastable phase density near a droplet will act on the intensity of the new 
droplets formation. This effect hasn't been taken into account in all previ- 
ous theoretical descriptions of the global evolution in the first order phase 
transition. But as it is shown in M this gap can lead to big numerical effects 
in the description of the whole process. 

The reason why the mentioned effect was not considered is rather trivial. 
Even under a spatial homogeneous consumption of the metastable phase 
the mentioned descriptions were rather complicate 0. Then ordinary the 
condensation process was described under the free molecular regime of growth 
where there is will be no such a gap. This was a serious restriction of the 
theoretical description. 

In some publications (see, for example f9|) the regime of the droplets 
growth was the diffusion one. This requires to consider the gap of the density 
near the growing droplet but in || the vapor consumption was regarded to 
be homogeneous in space. As far as this effect is very important one can 
not present a reliable description without taking it into account. Here we 
shall present a more realistic picture of a phase transition which allows an 
analytical solution. 

Qualitatively the picture of the condensation process is rather simple. 
A process of nucleation (i.e. a formation of supercritical embryos of liquid 
phase) leads to a vapor exhaustion which stops the nucleation process. But 
the supercritical embryos continue to consume a vapor phase. All surplus 
substance of a metastable phase will be accumulated in the embryos of a 
new phase. One can say that the process of condensation is completed^. 

The global picture of homogeneous condensation with explicit account of 
the density profiles was presented in |J where the giant numerical effects were 
observed. But ordinary the process of nucleation occurs on heterogeneous 

1 A further evolution includes the consumption of some relatively small embryos by 
some relatively big ones. It will be seen later that when all surplus substance is consumed 
then all droplets have approximately one and same size and we don't analyze this process 
here. 
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centers^. This radically complicates the theoretical description due to the 
centers exhaustion which has essentially non-linear character. 

During the nucleation process one can see that heterogeneous centers be- 
come the centers of the supercritical droplets which are growing irreversibly 
in time. But the nucleation process diminishes the number of free (unoccu- 
pied by droplets) heterogeneous centers. In some cases the total exhaustion 
of free heterogeneous centers interrupts the nucleation, in some cases the 
partial exhaustion seriously influences upon the nucleation rate. This effect 
has to be also taken into account. 

A simple analytical description of heterogeneous condensation will be 
presented here with proper account of all mentioned problems. As the result 
all main characteristics of the condensation process will be expressed through 
the parameters of external conditions behavior and substance parameters by 
some explicit analytical formulas. The error of the presented description will 
be estimated. 

The structure of the theory will be as follows: 

• At first we shall analyze the density profile around the solitary droplet 
and construct some approximations. It has much in common with 
the case of homogeneous condensation considered in || and will be 
considered briefly. 

• Then we shall construct some models for kinetics of process. We have 
to show that these models estimate the time evolution of the system 
during the nucleation period. As far as these models will give the 
similar results one can state that the approximate description of the 
nucleation kinetics is given. The error of description is, thus, estimated. 

• When the solution is already obtained we can compare it with formulas 
given by the previous approach and see the numerical effect of account 
of the gap near the growing droplets. 

The small parameter of the theory will be the inverse number of molecules 
inside the critical embryo of a new phase. The small value of this parameter 
isn't a restriction of our theory - it comes from the possibility of thermody- 
namic approach to calculate the free energy of the critical embryo. There is 
no other reliable way to calculate the free energy except the thermodynamic 

2 Also it is more simple to observe heterogeneous case experimentally. 
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approach^. To use the thermodynamic approach it is necessary to have at 
least a few dozens of molecules inside the embryo. 

Also we shall require the barrier character of nucleation. It means that 
to begin to grow irreversibly every embryo has to overcome the activation 
barrier of essential height. This height is less than for the homogeneously 
(pure fluctuation) formed embryo but still attains several thermal unitsQ. 
Certainly, one can imagine the situation when there is no activation barrier. 
Then all embryos begin immediately to grow irreversibly, The number of 
droplets (i.e. the irreversibly growing embryos) will be equal to the total 
number of centers and the kinetics of the process will be relatively simple. 

We shall speak only about the density profiles around the droplet and ig- 
nore the heat extraction in the nucleation process. Really, the mathematical 
structure of a diffusion equation is absolutely identical to the structure of a 
heat transfer equation. So, all constructions for the condensation heat ex- 
traction will be the same as for the substance consumption. It will lead only 
for some renormalizations. That's why only some corresponding remarks will 
be given. Any detailed results can be found in |10[ . 

We shall consider the situation of a metastable phase decay. It means 
that in the initial moment of time all substance is in the vapor phase. All 
heterogeneous centers are free from droplets. 



1 Profile around the solitary droplet 

Due to the external influence in the initial moment of time one can observe 
in the system the homogeneous mother metastable phase with the particle 
number density n which is equal to some value no- All heterogeneous centers 
are distributed rather homogeneously in space with the number density r] tot . 

The process of condensation can begin only when n Q is greater than the 
molecule number density in the vapor saturated over the plane liquid. 
The power of metastability of vapor will be characterized by the value of a 
supersaturation ( defined as 



3 A11 microscopic models requires very complex calculations which can not be fulfilled 
directly. 

4 All energy-like values will be measured in thermal units. 
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The initial value of the supersaturation will be marked by Co- 

Practically immediately around every center there will be formed an equi- 
librium embryo which has v e molecules of the condensing substance. The 
value of v e is relatively smallQ and there is no need to consider the density 
profile around the equilibrium embryos^. 

The number of the equilibrium embryos in the unit volumef] will be equal 
to i] tot . During the condensation process the number of the equilibrium em- 
bryos r\ will be decreased due to the exhaustion of the free heterogeneous 
centers 

V = Vtot ~ N 

where N is number of the supercritical embryos which will be called as the 
droplets. 

The last equation illustrates the main specific feature in kinetics of het- 
erogeneous condensation. Despite the simple form of the last relation the 
effect is very complicate because N depends on time in a very complicate 
manner. 

The effects of the density profile will be essential also for N and one can 
not directly apply the results of ||. One has to determine the effect of the 
influence on N even for the density profile of the solitary droplet. 

We shall call an approach where the law of embryos growth is found from 
continuous model but there is no account of the profile around the droplets 
as the "Additive approach" (AA). 

Then one can formulate the evident 

Statement 1 

The duration of the nucleation period^ and the characteristic sizes of the 
droplets at the end of nucleation period are greater than those calculated in 
AA 

5 In comparison with the characteristic number of molecules inside the droplet during 
the nucleation period. 

6 In fact it will disappear rather fast. This leads to the slight variation of the equilibrium 
embryo characteristics. This variation will act on the gap. But the final relaxation will be 
rather rapid. 

7 The system of a unit volume will be considered. 

8 The period of nucleation is the period of the relatively intensive formation of the 
droplets. It can be proven that the end of this period is well defined due to the cut-off of 
the intensity of the droplets formation. 
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Really, the existence of the density profile means that the part of sub- 
stance is going to be consumed from the regions where there is already no 
droplets formation. This substance is consumed from the gap instead of 
unexhausted regions as it is done in AA. 

Then having repeated all constructions^ from |J one can see0 

Statement 2 

The characteristic size of the droplets at the end of the nucleation period 
is many times greater than the size of the critical embryo. The main role in 
vapor consumption is played by the supercritical embryos. 

and 

Statement 3 

The characteristic time of the nucleation period duration is many times 
greater than the time of relaxation to the stationary state in the nearcritical 
region. Then one can use the stationary rate of nucleation as the intensity 
of the droplets formation in every current moment of time. 

Due to the statement 2 one has to investigate the profile around a growing 
droplet. The problem is whether one has to consider the interference of 
profiles around different droplets. 

To solve this problem one has to use the small parameter of the theory. 
Due to the statement 3 the rate of nucleation is equal to the stationary one. 
It can be taken from M 

I s = Zr)exp(-AF) 

where AF is the height of activation barrier (taken in thermal units), n is 
the number of free (unoccupied by the supercritical embryos) heterogeneous 
centers, Z is Zeldowitch factor. Zeldowitch factor is the smooth function of 
the supersaturation which is given by 

z - W 

ix l l 2 Av e Av c 

9 In [| the AA was formulated for external conditions of dynamic type. For the situation 
of decay the required hierarchical inequalities can be proven by the same way. One can 
note that in || there is no special reference on the type of condition when the required 
estimates are proven. 

10 The barrier character of nucleation is required here. It means that the magnitude of 
the activation barrier height has the same order as the free energy of the homogeneous 
critical embryos (may be it is 4 or 3 times smaller). 
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where W is kinetic factor, Av c is the halfwidht of the nearcritical region, Az/ e 
is the width of the equilibrium region. 

Due to the rather small size of the critical embryo it is reasonable to 
use the free molecular regime of the substance exchange^. Namely, in this 
regime an expression for the nucleation rate is well based. One has also 
to note that the critical embryo is in equilibrium (unstable one) with the 
metastable phase which implies no profiles of vapor density and regime of 
the substance exchange is the free molecular one. 

Under the free molecular regime W can be calculated as 

W = 3^u 2 J 3 a 

T 

where v c is the number of molecules inside the critical embryo, a is the 
condensation coefficient, 

is the characteristic time, v\ is a volume per one molecule in a liquid phase, 
vt is the mean thermal velocity of a molecule. 

The value of Av c is the halfwidth of the nearcritical region and it can be 
calculated as 

Au c = Yl exp(-F c + F^ir- 1 ' 2 

u<(u c +v e )/2 

where v is the number of the molecules inside the embryo, F v is the free 
energy of the embryo of v molecules, F c is the free energy of the critical 
embryo. In continuous approximation it can be calculated asQ 

Av = I — I l 1/2 

The value of Az/ e can be calculated as 

Av e = £ exp{-F v + F e ) 

u<(u e +u c )/2 

11 As long as the characteristic size of the droplet during the nucleation is many times 
greater that the critical size the diffusion regime of growth for the characteristic droplets 
is quite reasonable. 

12 Ordinary Av e is smaller than Ai/ C and an explicit summation for lS.v c is quite 
reasonable. 
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where F e is the free energy of the equilibrium embryo. 

Both Au c and Az/ e are rather smooth functions of the supersaturation. 

One can see that I s is a very sharp function of the supersaturation. It 
means that the relatively small fall of the supersaturation leads to the inter- 
ruption of the droplets formation. 

At least for all ( > Co/2 one can show that d 2 (/dt 2 > and the is no long 
tail of a size spectrum with a small intensity. It means that the interruption 
of the droplets formation leads to the interruption of the nucleation process. 
So, the relative fall of the supersaturation during the nucleation process is 
small. One can come to the following statement 

Statement 4 

During the nucleation period the relative variation of supersaturation is 
small. 

The last statement shows that there is no need to consider the interference 
of profiles in order to change the rate of droplets growth. 

On the base of mentioned expressions and the smallness of the relative 
fall of supersaturation one can see the validity of the following approximation 

J,(C) = / s (Co)exp(AF(Co)-AF(C)) 

for the nucleation period. 

Moreover one can linearize the height of the activation barrier over the 
supersaturation and get 

/.(C) = /.(Co) ex P (-^% =Co (C - Co)) (1) 

The validity of the last approximation depends on the concrete type of het- 
erogeneous centers but for the majority of the heterogeneous centers types 
the last approximation is valid. For example, this validity can be directly 
proven for ions. 

One can explicitly calculate the derivative in the last expression^ 

dAF 1 , N 

—77- = --r—{yc-v e ) 
d( C + 1 



13 Here we assume the vapor to be an ideal gas and suppose the possibility to present 
the free energy of critical and equilibrium embryos as an analytical function of an inverse 
embryo radius. 
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The smooth character of the last expression shows the validity of (|1|) one 
more time. 

Then ([I]) can be rewritten as 



/ s (C) = / s (Co)exp(r 



C-Co 



(2) 



Co 



where 



r = -Co 



dAF 



Co 



^e(Co)) 



C=Co - 



Co + 1 



As far as the value of v c in going in the presented theory to infinity the 
value of T is also going to oo. The real value of T is very big. Certainly, 
one can consider the possibility of compensation between v c and v e in the 
expression for T. Then one has to mention thatdue to the barrier character of 
nucleation at least v c — v e > Av c . Having estimated Av c by the homogeneous 
value Av c ~ z/?/ 3 one can see that T 3> 1 in any case. 

The small value of T" 1 will be very important in further constructions. 

We see that essential dependence over the supersaturation occurs through 
the height of activation barrier. It allows to give the interpretation of the sta- 
tionary rate of nucleation as the probability for the given embryo to overcome 
the activation barrier. After the interpretation of l s as the probability we can 
apply it to the arbitrary spatial point of spatially unhomogeneous system. 
To use this interpretation the natural requirement is the weak unhomogenity 
of the system. Namely, the volume of the regions where 



is violated has to be relatively small. 

Here D is the diffusion coefficient, t s is the time of relaxation in the 
nearcritical region which can be estimated according to Zeldowitch |l| 



One can use instead of t s the time Z which can be interpreted as the mean 
time to overcome the nearcritical region. 



C(r) - C(r + y/4DQ 
C(0 



-i 
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Both these estimates are observed. Actually we need them only for those 
regions where the intensity of the droplets formation isn't too small in com- 
parison with initial intensity. Certainly, the required property is observed 
there. 

Now we have to turn to the rate of embryos growth. According to the 
Statement 2 the characteristic size of droplets is rather big. Then it is more 
reasonable to use the diffusion regime of droplets growth. Under the inter- 
mediate Knudsen numbers one has to use an interpolation law for the rate of 
embryos growth (for example, see [|J). It will be important that all expres- 
sions for the embryos growth lead to the avalanche character of the substance 
consumption. 

The avalanche character of the substance consumption means that the 
quantity of substance accumulated by a droplet strongly increases in time. 
The most evident manifestation of the avalanche consumption can be seen in 
the free molecular regime of the substance consumption. The most "weak" 
effect of the avalanche consumption can be seen under the diffusion regime 
of the substance consumption. 

The force of the iterations convergence in || is based on this property. 
The property of avalanche consumption will be extremely important in the 
further constructions also. That's why we take the diffusion regime to have 
the worst situation and to grasp errors in all possible cases. 

In the diffusion regime of the vapor consumption the law of growth for 
a droplet (i.e. for the supercritical embryo) can be written in the following 
way 

dt s 



where 

K 



is some constant. 

The last expression is written in stationary approximation. The non- 
stationary effects were investigated in many publications in details and they 
are rather small ones. 

One can see that the rate of the droplet growth is proportional to £. So, 
the rate of growth can be essentially changed only by the essential relative 
variation of (. Then according to Statement 4 one can see that 

Statement 5 
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The rate of droplets growth during the nucleation period can be approxi- 
mated as a constant one. 

The last statement is extremely important because is allows to analyze 
the profile of the density initiated by a solitary droplet. Now this case will 
be considered. 

The approximately constant value of the supersaturation allows to inte- 
grate the law of growth and get 

u(t) = 7 * 3/2 

where 

7 = (4 vr )3/ 2 (^) 1 /2(^22^)3/2 

and t is the duration of the irreversible growth for the given droplet. 

Consider the spherical system of coordinates with a center in the center 
of droplet. Diffusion equation will be written as 

dn ^ . 
— = DAn 
dt 

where A is the Laplace operator. Diffusion coefficient D is supposed to be 
approximately constant (there is a lot of a passive gas and the density of a 
gas mixture is approximately constant). 

The boundary conditions are the following 

n\ r=OQ = ra(oo) 

n\ r =R d = Uoo 

where Rd is the radius of a droplet. The values and n(oo) are known 
parameters. The variable r is the distance until the center of the embryo. 

The stationary approximation is suitable for the rate of droplets growth. 
The errors are analyzed in |4j] and they are small. But the stationary solution 
can not give any reasonable result for the density far from the droplet. Really, 
the stationary solution is the following 

R 

n{r) = n(oo) - — (n(oo) - n^) (3) 
and has a very long tail. This tail leads to the infinite value of 

poo 

G = 47rr 2 (n(oo) — n(r))dr 
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which has to be the integral excess of the substance. This quantity has to be 
in the droplet. This contradiction shows that it is absolutely impossible to 
use the stationary approximation for the density profile around the droplet. 
One has to introduce another approach. 

One can see that if the first boundary condition will be changed by 

n\ r=OQ = n(oo)(l - r _1 ) 

then the rate of embryos growth will be changed unessentially. But the level 
n(oo)(l — r _1 ) is the level when the nucleation stops. So, one can see that 
during the nucleation period there is no interaction between droplets through 
the change of the growth rate. Certainly, two droplets can appear too close 
and act one upon another but the probability of such coincidence is too small. 
That is why one can come to the "Principle of a separate growth of droplets 
during the nucleation period" . 

Now one has to prove that at the distances (5 -f- 10) Rd from the droplet 
one can observe the quasistationary profile. 

One has to note that 

v t /v v < 1 (4) 

where v v is the partial molecule volume in a vapor phase. Really, the last 
ratio is very small (for example, this ratio is 0.001 for water in normal ther- 
modynamic conditions). But contrary to T" 1 one can not consider it in all 
cases as zero. 

Now one can introduce a formal parameter I which attains big values 

/ > 1 

but satisfies the following condition 

l 2 ^ « 1 (5) 

v v 

Due to (£|) it is possible to do. 

In the region r < lR d the stationary profile is established after 

_ l2 Ri 

It is necessary to show that 

R d (t + t h ) - R d (t) 

s " — m — <<x 
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Really, 

dRd th 
s ~ 



dt Rd 
and 

,2 V l 

s ~ I — 

v v 

which is a small value according to @. So, the stationary form of profile in 
the region r < Rdl is proven. 

As far as T ^> 1 and at least V ^> / one can see that in the region r < IRd 
there is no formation of new droplets. Then this region isn't interesting for 
the theory and one can only observe the region r > IRd- 

The previous notation is rather important. Namely this property allows to 
use the model with a point source. Really, one can consider only the distances 
greater than IRd- But at these distances the droplet can be interpreted as a 
point. 

Certainly, the point approximation of a droplet can not give an expression 
for the rate of droplets growth because the boundary condition at r = Rd is 
absent. But the rate of growth is already known and can be used directly as 
a known function of time. Really, 

& = A *V2 
dt 

where 

The action of a point source of a vapor consumption can be described 
in a simple and suitable manner under the Green function formalism. The 
Green function for the diffusion equation can be written in a following form 



Gr = G(t) 



(47rL>t) 3 / 2 

Then one can get the density profile by a simple integration 



2 



f \ ( \ f* XX1/2 ( T ■ / 

n(r) — n(oo) — / — — — — — exp — — — -)dx 
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After evident transformations one can come to 



where 




and 

/(£)= / (^--) 1/2 exp(-x 2 )dx 

Jj3 (3 Z X z 

It is important that the profile dependence on t and r is now going through 
(3- 

Concrete form of f(0) is drawn in fig. 1. One can get for f{(3) an 
expression through special functions 

/(/3) = ir(^)ex P (-m!^/? 2 ) 

Here V is the Gamma-function , \I/ is the confluent hypergeometric function. 

One can get asymptotes for f(/3) at small and big values of (3. At small 
values 

m ~ % m 

which corresponds to the stationary solution fl3|). 
At big values of (3 one can come to 

f{(3) = exp(-/? 2 )^^ °x 1 / 2 exp(-x)^ ~ (8) 

One can see that this asymptote radically differs from the stationary solu- 
tion. Namely this tail behavior ensures the convergence of the integral for G. 
Certainly, the Green function formalism ensures a precise value for G which 
is introduced here as an external object. 

Now one can turn to construct some approximation for the nucleation 
rate around the growing droplet. 

One can see that according to (0) the behavior of the supersaturation is 
important when (q — ( < (2^3)Co/T- When ( — ( > (2^3)( /T the intensity 
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of the droplets formation is negligibly small. From ([!]) one can see that 



J s (C(r))=/4Co)exp(-rJ^/^/(/3)) 




Then one can extract a positive parameter 




which will be important in further constructions. 



Due to T ^> 1 one can easily see that 



a > 1 



The last condition isn't necessary for further constructions, but it will be 
rather important for manifestation of profile effects in the nucleation process. 

The last condition is also the most doubtful one because v\jv v <C 1 and 
one has the ration of two big parameters with generally unknown result. It 
is necessary to stress that condition v\jv v <C 1 isn't so strong as T 3> 1. In 
frames of thermodynamic description T 3> 1 is the main condition necessary 
to give the thermodynamic description and v\jv v <C 1 is the supplementary 
condition which isn't absolutely necessary but slightly simplifies the theory. 

In the situation of homogeneous condensation one has the "hidden con- 
tradiction between the thermodynamic description and the relatively in- 
tensive nucleation". Really, as far as in the homogeneous condensation 
AF = F c ~ v^ 3 the limit v — > oo means AF — > oo and the rate of nu- 
cleation goes to zero. So, there is the contradiction between thermodynamic 
limit in the critical embryo description and the observable^ rate of nucle- 
ation. 

In the case of heterogeneous condensation there is no such contradiction 
if there are some active centers of condensation. Then the height of the ac- 
tivation barrier has no direct connection with a number of molecules inside 
the critical embryo. For example, the halfwidth of the nearcritical region es- 
timated by a homogeneous value has a value ~ z/ 2//3 and goes to infinity when 
v — ► oo, but the free energy decreases at the boundary of the nearcritical re- 
gion only by one thermal unit. So, in certain sense the case of heterogeneous 
condensation is more preferable for theoretical description. 

14 Which isn't too small. 
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As a compensation for this advantage one has to note that the Statement 
1 and Statement 2 are based on the homogeneous estimate for the activation 
barrier height. These properties can be violated. But as far as these State- 
ments are based on very strong inequalities one can accept their validity. The 
influence of nonstationary effects will be analyzed in a separate publication 
where the complete theory has been constructed. 

Now one can analyze the profile of intensity of the droplets formation 
around the already existing droplet. This nucleation rate profile is a rather 
sharp function which has a step-like behavior. 

To show this we shall introduce two characteristic values of j3 (/3 s t and 
(3 fin) by relations 



In the region (5 > (3 s t the rate of nucleation practically coincides with the 
unperturbed value I s {Co)- In the region ( < Q in the rate of nucleation is 
negligible in comparison with the unperturbed value, i.e. I s (((r)) <C I s (Co)- 

In some moment t the values f3 st , (3f in initiate the space distances r st , rf in 
by the following expressions 





P fin V4JDt 



For (7> 1 one can come to 



/(£*) « i 



f (13 fin) < 1 




fist + P fin ^(3stjin 
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\ r st-r fin \ = 1 ^ 1 
r s t + r fin A(3 st)} in 

The real picture of nucleation is going in the time scale. In a fixed space 
point r one can introduce two characteristic times t st and tf in by expressions 

r 2 

Before t st one can not observe any essential deviation of the nucleation 
rate from the unperturbed value. After tfi n the rate of nucleation is very 
small. 

One can get for the relative deviation 

^ tfin ^st 



the following expression 



tst, fin 



1 



Pst,fin 

So, the relative deviation is small. 

Even in the situation of small a one can show with the help of the asymp- 
tote (0) that 5 is rather small. 

The step-like behavior^ of intensity profile allows to introduce some char- 
acteristic parameter (3 e ff and consider the regionF] 

P < Peff 

as the exhausted region where there is no nucleation more and in the region 

P > Peff 

the rate of nucleation is unperturbed^. 

15 This step can take place even for /3 > j3 s t when a -C 1. 
16 At a <C 1 the value of /3 e // can be greater tan fist and /3/i„. 
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In all cases (3 e ff > f3fi 



17 



One has to choose (3 e ff carefully. The problem is the possibility of ex- 
istence of a long tail of a density profile. To grasp the situation of a small 
values of a one has to introduce f3 e ff in an integral manner. 

One can introduce the excess of nucleation rate AI S by the following 
formula 



when I s is the unperturbed rate of nucleation. 

On the base of the last expression one can get the excess of N due to the 
existence of the solitary profile. This value will be marked by AN so i and can 
be found as 





Having used equation (H) one can come to 





The value AN so i can be now found as 



AN sol = A7c(4Dt) 3/2 L 



rt poo 

s / (l-exp(-r 
Jo Jo 




f(p)))p 2 dpdt' 



The step-like approximation of the nucleation profile will lead to 



Jo 



The value (3 e ff has to be determined from 



{(3 e}f ) = AI S 
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Inside the volume 

4 n 

Veff = gvrr e// 

there is no nucleation and outside this volume the rate of nucleation is un- 
perturbed. 

Then one can imagine that around every droplet there exists the ex- 
hausted region (ER) where there is no nucleation and the unexhausted region 
(UR) where the nucleation remains unperturbed. The whole space now is 
divided into two regions. 

The volume V e ff grows in time in the following way 

Veff = fn(3t ff D^ 

In a free molecular regime V e ff will grow even faster. 
For b e ff one can get the simple expression 

fi ff = 3 f°(l - exp(-r JI/^/(/?)))/? 2 d/? 

13 Jo V V v v 

or 

P! ff = 3^°°(l - eM-v 1/2 ^Jlf(P)))P 2 dP 
For AN so i one can obtain 

AN sol = I s (( ) f dt'V eff = fdt'Ur 3 ^ 
Jo J o o 

One can easy integrate the last expression and get 

One can see that AN so i is growing in time rather rapidly. Namely this 
property illustrates the feature of the "avalanche consumption during the 
first order phase transition" in application to the heterogeneous nucleation. 
For those situations where tr>l one can get 

Peff ~ Pst ~ /3 fi n 
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and p e ff satisfies the simple equation 



The last equation can be easily solved by iterations as far as /3 e // 3> 1 and 
exp(— f3 2 ) is the most sharp function of the supersaturation. 

Earlier when the principle of the separate growth was discussed then the 
reason of the absence of interaction between droplets was the low probability 
to appear too close one to another only due to the smallness of such space 
volume. Now one can see that the growing ER helps to exclude the inter- 
ference. Really, the essential deviation of the supersaturation from the ideal 
one can be seen in the region r < R^l. It means that the distance between 
the droplets with an interference must have the order 2RJ, Then the time 
distance between the moments of formation of these droplets must be shorter 
than 

At mtt ~ (^t^) 

This time interval is many times shorter than the duration of the nucleation 
period. 

Every droplet forms rather rapidly after formation the ER of such a size 
which guarantees that the rate of growth of the given droplet can not be 
perturbed by a vapor consumption by the other droplets. 



2 Kinetic models of the global evolution 

Now one can construct the picture of nucleation in a whole system. The 
main problem is to take into account the interference of the density profiles. 
The interference through the rate of growth is absent, but there is a simple 
overlapping of profiles which leads to the deviation of the total nucleation rate 
over the volume from those calculated with account of the additive excesses 
around every droplet. 

The overlapping of ER (even when this approximate formalism is used) is 
very complex and can not be directly taken into account in precise manner. 
Instead of writing the long expressions which can not be explicitly calculated 
one can act in another manner: some simple approximate models for kinetics 
of the nucleation process will be formulated. From these models it will be 
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seen that they estimate the nucleation characteristics from below and from 
above and lead practically to the similar results. So, it will be shown that 
the complex details of the ER overlapping has no strong influence on the real 
characteristics of the phase transition. 

At first one can consider the common feature of all models concerning 
the exhaustion of free heterogeneous centers. 

The rate of nucleation I depends on time t and on spatial point r (the 
last behavior is the most complex). So, it is reasonable to consider the mean 
(over the space) value of I denoting it by < J >. For < I > one can write 
the following expression 

< / >= ^=JL, (9) 
W tot r] tot 

where Jo is the unperturbed rate of nucleation. Here Wf ree is the volume of a 
region where the rate of nucleation is unperturbed, i.e. the UR of the whole 
system. The value Wtot is the total volume of the system (it equals to unity 
and is written only to clarify the consideration). 
Then as far as 

N = [ <I> (t')dt' 
Jo 



one can get 



r < / > (t')dt' 

Jo 



V = Vtot 

In the differential form the last relation can be written as 

dr] 

— L = -<I> 
dt 

and with the help of (|9|) it can be rewritten as 

dV _ V Wfree. j 

dt Vtot W tot ° 
After the integration of the last expression one can come to 



7] = ?7to t exp( 



Wtot Vtot 



One has to note that heterogeneous centers aren't distributed homoge- 
neously with respect to the ER (or UR). Only free heterogeneous centers are 
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distributed homogeneously with respect to ER. This fact has to be also taken 
into account. 

The problem is to determine the value of Wf ree . In different models it 
will be determined in different forms. 
The model without overlapping. 
One can write 

w free = w tot - w exh 

where W ex h is the volume where there is no further formation of the droplets. 
Very approximately one can present it as the sum of ER around all appeared 
droplets 

W exh « J2 V eff 

i 

(the sum is taken over all already formed droplets). Certainly, the last ap- 
proximation is rigorous only when there is no overlapping of the ER around 
different droplets. 

Having used an expression for V e ff one can come to 

W exh = J*dt'<I> (t')f n(3 3 eff D 3 /\t - t'f' 2 (10) 

Having used an expression for < / > one can come to the closed system 
of nucleation kinetics equations 

W free = W tot - fdt'^^I*X(3l ff D^{t - t'f' 2 
Jo Vtot Wtot 3 

77 = r] tot exp - / — j- dt 11 

Jo W tot Vtot 

In the quasihomogeneous limit (when there is no essential exhaustion of 
centers) this system can be reduced to 

W free = Wtot - f dt'^I*Xpl n D^{t - t'f 12 
Jo Wtot 3 

which can be written after the evident renormalization t — > at, t' — > at' where 
a = (I Y^PeffD 3/2 ) 2/5 in the universal form 



W free = 1 - / dt'W free (t - t') 

J 



A 3/2 
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One has to note that in the general case the system of nucleation equations 
can be solved with the help of methods developed in 0. At first one can 
solve the quasihomogeneous equation (it the Volterra equation with a rather 
simple kernel which allows to apply the Laplace transformation to solve it) 
and then on the base of quasihomogeneous equation one can find the final 
rather precise expression using the second equation as a formula for r\. 

Another variant is to solve numerically the universal equation for Wf ree } lom \ 

Wf ree horn = 1 ~ J dtfW free hom(t ~ t') 3 ^ 2 

Then one has the universal function Wf ree n om- Then one can find rj as 

. i /"* Wf ree homit) Iq ,,/s 

7] = rjtot exp(-a / ... at ) 



10 Wtot Vtot 

The last expression leads to the formula for < / > as: 

< / >= WfTee hom expt-a- 1 I* Wfree —dt%) 
Wtot ' Jo W tot Vtot 

The justification of such approach is analogous to 0. The physical reason is 
very simple: when there is no exhaustion of heterogeneous centers then the 
solution is found precisely, when there is an essential exhaustion of centers 
there is no need to know Wf ree with a high precision because the converging 
force of the operator of the right hand of the second equation in the nucleation 
equations system is extremely high. 

Now we shall take into account the overlapping. It can be done rather 
approximately. 

The model with chaotic overlapping 

The matter of discussion is the correct expression for Wf ree which can not 
be got absolutely precisely. Now the more reasonable expression for Wf ree 
will be presented. Certainly this will lead to a more complicate equation 
which will be more difficult to solve. 

One can use the differential approach to write the expression for Wf ree . 
Having written the evident relation 

dWfree _ dW exh 

dt ~ dt 
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one has to invent an approximation for dW ex h/dt. Here the following approx- 
imation 

dW exh _ dJ2jVeffW free 

dt ~ dt Wua 

(the sum is taken over all droplets) will be used. It corresponds to the 
following approach: The probability of the absence of overlapping of the new 
parts of ER around the given droplet with other ER is proportional to the free 
volume of the system. 

The last supposition seems to be rather reasonable. 

The value dJ2i Veff/dt can be rewritten as 

d T.jVe}f _ dVeff 

dt i dt 

The last value can be easily expressed through < I > as 

X^$ L = \f Q * f<I> Mf*#//tf 8/2 (* " 1/a (12) 



due to (0). 
Then 



<IU '" h 3 ft < I > (t')^ ff D^(t-t')^dt'W free (t) 



IT' 2 Jo ^ ^,y^ e// 

and 

dW free 3 /•' . ,32 



— - -- Jo < I > (Oy^e//£> 3/2 (* - t') 1/2 dt'W free (t) 
Having used an expression for < I > one can come to 



Together with ( |TT] ) the last equation forms the closed system of the nu- 
cleation equations in the second model. 

The previous equation can be integrated which gives 

In W free = -f Eg££^l /o I^ 32 £ D 3/ 2(t _ £ / )3 /2^ + CQnst 
Jo W tot r) tot 3 JJ 
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Due to the initial conditions the const is equal to zero. Having intro- 
duced the function F = — In Wf ree one can get for F the following system of 
equations 

F(t) = r eX p(-F(t'))/o^f ^l n D^{t-t'f /2 dt' 

JO t]tot o 

ft J Q 

7] = 7] tot exp(- exp(-F(t)) — dt) 

JO T] t ot 

One can see that the system of condensation equations is identical to 
the system of nucleation equations in AA. It was completely analyzed in |J. 
Certainly, the parameters in the system will be another. 

The last system can be rewritten after the evident renormalization as 

F(t) = I* exp(-F(lf))(t - t'f /2 6(t')dt' = F(F, 9) 
Jo 

6{t) = exp[-A f exp{-F(t'))dt'] = 6(F) 
Jo 

where 8(t) = rj(t)/rj to t and A is some known parameter. 
This system can be solved by iterations defined as 

Fi+i = F(Fi, 6i) 

Qi+i = Q{Fi) 

with F = 0,9 = 1. 

For Fi, 9i one can get the chains of inequalities 

F < F 2 .. < F 2l < ... < F < ... < F 2i+1 < ... <F 3 < F 1 

e 1 <9 3 < ... < e 2i+ i < ... < e < ... < e 2l < ... <e 2 <e 

Then one can estimate errors in F^, 6i. 

One can use also another methods analogous to those described in [[J. 

The similarity of the the condensation equations in AA and in the second 
model is extremely important for the transition towards collective character 
of vapor consumption which is analyzed in ||. 

The physical ground of the considered model is the chaotic overlapping 
of ER. Namely, the chaotic overlapping lies in the base of approximation 
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used here. But due to the spherical form of every ER the overlapping isn't 
absolutely chaotic. 

What can be done in such a situation? In the next model we shall show 
that concrete type of overlapping isn't very important. 

To finish with the second model we shall show the same method of its 
solution as that of the first model. 

One can also formulate the quasihomogeneous equation as the following 

one 

r t 32 
F hom (t) = j exp(-F hom (t'))I —Trf3 3 eff D 3 / 2 (t-t'f 2 dt' 

Then rj can be approximately found as 



V = Vtot 



exp(- / exp(-F hom (t')) dt 

Jo r] tat 



The quasihomogeneous equation can be renormalized. After renormal- 
ization z —>■ at, t' — > at' where 

a = (/of ^£3/2)2/5 

one can transform the quasihomogeneous equation into the universal form 
\nW free hom {t) = - f W free hom {t'){t - t'f /2 dt' 

^ 

The model with the formation of droplets inside the ER 

The third model will show that the role of overlapping isn't so essential 

as it can be imagined from the first point of view. 

Suppose that the new droplets can appear also in the ER of the already 

existed droplets. Then instead of (O) one has to use 



2 Jo rjtot 3 

Then 



dW exh 
dt 2 Jo rjtot 



3 r* m io ^ ffD 3/2 {t _ t iy/2 dt i Wfree{t) 

JO That O 
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and 

dW free 3 ft 77(f) 32 3 3/2 ,/xl/2,,/ w /.n 

^T" = — 2 J Q h — -^ en D ' (t — t) ' dt W free (t) 

Together with (JTTj) the last equation forms the closed system of the nucleation 
equations in the third model. 

The substance balance equation of the system can be integrated which 
gives 

In W free = -( Io^l-Tr^ffD^it - t'f' 2 dt' + const 
Jo ritot 3 JJ 

Due to the initial conditions the const is equal to zero. Having intro- 
duced the function F = — In Wf ree one can get for F the following system of 
equations 

F(t) = fl Q ^X{3l ff D^{t-t'f ,2 dt' 

JO Vtot "J 



Vtot 

7] = r] tot exp(- I cxp 



F(t'))—dt') 

Vtot 

This system corresponds to the first iteration in the iteration solution of 
the second model by the method of iterations of a type "a" described in [|]. 
It is also described here. But for rj the whole set of iterations has been taken 
(see details in ||). 

One can slightly modify the model and suppose that in the expression for 
n one can use the same approximation for < / > as in equation for Wf ree . 
Then the system of equation will precisely correspond to the first iteration 
in the iteration solution in 0. 

Actually one needn't to analyze these models in details following || but 
has only to note that all these solutions are practically similar. One can see 
the solution of the third model (both variants) after the same renormaliza- 
tions as in the second model. 

Now one has to explain why the third model is rather accurate. One can 



do it only with the help of results obtained in JLO . 

In JTIJ it was noted that when the power of the kernel (t — t') is rather 
big then the solution of the quasihomogeneous equation weakly depends on 
concrete value of a power. 

Another important notation shows that when the power of (t — t') is 
extremely high then the ER of the given (first) droplet formed inside the 
ER of another (second) droplet practically can not go outside the ER of that 
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(second) droplet. Then the third model is absolute adequate in this situation. 
The same feature can be seen directly from the iteration procedure results. 

Having combined these two notations one can see that on one hand the 
second model is close to the situation with great power of the kernel and in 
the situation with the big kernel the third model is suitable. 

Now it is possible to explain why the overlapping isn't so important as 
it can be imagined. Really, as far as the power of kernel is big and one 
can observe the avalanche consumption one can see the following qualitative 
picture 

• Practically during all period of nucleation the total ER is small and 
there is no problem of overlapping. 

• At the end of the nucleation period the total ER is going to occupy 
the essential part of the volume and a few moments late it occupies 
all the volume of the system. This process is very rapid. It stops the 
nucleation. 

This picture shows that there is no strong influence of the overlapping on 
the nucleation process (except the cut-off). 

As the result the nucleation description is completed. One can use both 
the second and the third model to get the nucleation description. How to 
solve these equations is also described here. 

The method of universal solution 

The main idea of the theory presented in || was to consider the quasi- 
homogeneous equation, to get the universal solution and then on the base 
of this solution to calculate the number of free heterogeneous centers. As 
the result one can get an expression for < I > and can calculate the total 
number of droplets appeared on the nucleation process. 

Here one can follow the same idea. But one can make this idea more 
deep. One needn't to formulate the universal quasihomogeneous equation. 
Instead of equation one can formulate the universal model. 

The model will be the following 

• -the rate of nucleation < I > can be found from 

< / >= h Ztp± JL 

Wfot Vtot 
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• - with intensity Iq the droplet appears in the arbitrary point of the 
system which belongs to the free volume. 

• -the value Wf ree can be found by exclusion of all ER around the already 
existing droplets. 

• - if the point is occupied by the ER of any droplet then the droplet can 
not be formed. 

• - the size r e ff of ER grows in time according to 

r eff = 2(3 eff D^ 

• -the initial condition is the absence of droplets and random distribution 
of centers. 

With the proper renormalization of time t and size r one can cancel all 
coefficients. Then this process will be the universal one and as the result 
Wf ree is the universal function of time. Then one can directly apply ([LI]) and 
get the number of free heterogeneous centers (after the proper renormaliza- 
tion)Q. 

The modification for the dynamic conditions requires to put instead of Iq 
the value Ioexp(ct) with some parameter c determined by external conditions 
and to change the lower limit of integration to — oo ||. 

The main constructions of the theory will be absolutely the same but the 
forms of characteristic curves changes radically. 

3 Numerical results 

Numerical simulation in the nucleation kinetics plays at least two important 
roles. The first role is the standard comparison with the approximate models 
to observe their quality. The second one is more specific and it is concerned 
with some universal dependencies in the nucleation kinetics. 

In the AA to the nucleation kinetics it was shown that the adequate 
approach can be presented on the base of the quasihomogeneous solution 
0. Despite the dynamic conditions considered in one can observe this 

18 This is quasihomogeneous approach method 
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property in the situation of the metastable phase decay also. Recall the 
reasons for such approach. 

The formal reason is the careful analysis of the iteration procedure pro- 
posed in |U . Really, the final result for the total number of droplets appeared 
in the nucleation process is given by the second iteration (see iterations of 
type "a" in ||) for the relative number of free heterogeneous centers. This 
iteration is based only on the first iteration for the supersaturation. The 
value of supersaturation there is calculated without any account of the het- 
erogeneous centers exhaustion. So, one can see that the final result can be 
obtained on the base of the supersaturation in the quasihomogeneous ap- 
proximation. This approximation can be more sophisticated than the first 
iteration. Namely this approach was used in []5| where the precise quasiho- 
mogeneous universal solution was chosen as base for the final results. 

The physical reason of such behavior is rather simple. The main role 
in the vapor consumption belongs to the droplets of relatively large sizes. 
We have already marked this fact. But moreover, due to the avalanche 
character of the vapor consumption the main role in vapor consumption is 
played by the relatively large droplets which are formed in the first moments 
of time of the nucleation period. When the effect of the centers exhaustion is 
essential already in the first moments of nucleation period^ then at the end 
of nucleation period already all centers will be exhausted. Then the result 
is evident - all centers are the centers of droplets. Due to the high force 
of convergence in this situation this result can be gotten under the rather 
arbitrary behavior of supersaturation (even including the quasihomogeneous 
case). In the opposite case the exhaustion of heterogeneous centers during 
the first moments of the nucleation period isn't essential and one can get the 
quasihomogeneous behavior of supersaturation. 

This effect will be called as the "approximate separation of the hetero- 
geneous and homogeneous problems". It is based only on the avalanche 
consumption of the metastable phase. So, there are no objections to see this 
effect also in the situation with the density profiles. Then it is extremely im- 
portant to get solution in the quasihomogeneous situation and clarify whether 
it can be presented in the universal form. 

19 More precisely one can define these "first moments of time" as 2/5 of the nucleation 
period duration (under the free molecular regime it is 1 /4 of the nucleation period dura- 
tion). The reason for concrete value can come from the iteration procedure. 
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The universal form of the quasihomogeneous solution can be easily seen 
in the situation of profiles also. In the AA there was no specific space scale 
because the consumption took place homogeneously from all space points of 
the system. Here in the situation of profiles there is the elementary space 
scale and one can choose the space scale as to have that the linear size of ER 
around the droplet is growing^] as t 1 ^ 2 . The time scale has to be chosen as 
to get that in the free volume equal to the total volume of the system one 
can see the appearance of one droplet in the unit of time. As long as the 
functional dependencies of the nucleation rate and the radius of ER are not 
identical one can do such a renormalization without any problems. So, we 
see that the pseudohomogeneous case allows the universal description. 

The process of the heterogeneous centers exhaustion destroys this univer- 
sality and one has to act as inQ ||. The total number of droplets has to be 
approximately calculated as 

Ntotai = Vtat[l - exp( — )] (13) 

Vtot 

where Nh om is the number of the droplets appeared in the quasihomogeneous 
situation (under the same parameters). This formula can be also used for all 
approximate models described earlier. 

For the numerical simulation it was convenient to consider the cube box 
with a size 10. The rate of the ER growth is chosen as 

dR i /o 

= 100t 1/2 

dt 

where R is the radius of ER. The rate of nucleation is chosen as to have one 
attempt of a new droplet formation in the system during dt = 0.002. This 
attempt is governed by the random procedure. It may lead to the point in 
one of the ER and then no droplet will be formed. In the opposite situation 
when the point is indicated out of all ER of the already existing droplets 
there will be a new formation of a droplet. 

One has to stress that the random procedure has one specific negative 
feature. In the standard numerical procedures the next random coordinate 
is calculated on the base of the previous ones. So, if the current coordinate 



20 Certainly, the power has to be conserved. 

21 Here we use slightly more simple way. 
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lies near the center of the already existing ER then the next coordinate 
will be also near the center of another ER. These correlations lead to the 
necessity to consider the large system. In the system under consideration 
the number of droplets appeared in the quasihomogeneous situation will be 
near 500. Nevertheless the square mean fluctuation will be about 20. The 
careful consideration can show that the error induced by the substitution of 
the periodic boundary conditions by the zero ones has the same power as 
the mean square error. It can be directly seen by the numerical simulation. 
It is explained by the evident notation that the characteristic overlapping 
of profiles is about the mean profile size. We shall call this feature as the 
property of "moderate overlapping". This fact can be proven analytically. 
One can also analytically prove that the effect of fluctuations is really caused 
by the mentioned phenomena. 

The mean value of the total droplets number is equal to 504.8 (under 
the zero boundary conditions). This value has to be put into the previous 
formula. 

The avalanche character of the vapor consumption is illustrated by fig. 
2-4. Three different moments of time t = 0.5, t — 1 and t = 1.5 are chosen 
as characteristic values. The cross section is drawn. The dashed regions 
correspond to the ER of the already existing droplets. The black regions 
correspond to the overlapping^ of ER. 

The number of heterogeneous centers in such a system can be arbitrary. 
Certainly the effects of the centers exhaustion will be important when the 
number of centers is small (in comparison with Nh om = 504.8). The pictures 
for r] tot = 50 are drawn in fig. 5-7 for t = 1.5, t — 3, t = 6. One can see that 
the number of ER is smaller than in the quasihomogeneous case. The size of 
ER when the free volume is going to be exhausted is larger than the size of 
ER in the quasihomogeneous case. The time necessary to cover all volume 
by ER is greater than in the quasihomogeneous case. It doesn't mean that 
the duration of the nucleation period will be longer (simply all centers will 
be exhausted and this means the end of nucleation). Moreover, the duration 
of the nucleation period in the situation of the relatively small number of 
heterogeneous centers will be shorter than in the quasihomogeneous case. 

One can also see that the avalanche character of the vapor consumption 
here will be more smooth than in the quasihomogeneous case. Certainly, in 

22 Only when the distance between centers is approximately the odd one. 
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the quasihomogeneous case the appearance of the new ER helps to consume 
the vapor phase in the avalanche manner. But in the situation with the small 
number of centers there is no need to consider the process carefully because 
the exhaustion of centers leads to the evident result of condensation - the 
number of the droplets equals to the number of centers. 

Now it is evident that the main object of our interest will be the quasi- 
homogeneous case. The relative rate of nucleation in this case is shown in 
fig. 8. Here the rate of nucleation is averaged over lOOdt = 0.2 and over 16 
attempts. So, the rate of nucleation is rather smooth function. 

The relative rate of nucleation is compared in fig. 8 with the mentioned 
models. Certainly, the rate of nucleation defines the spectrum of sizes when 
we take as the size of embryo those characteristic which has the rate of growth 
independent on the value of size. For the diffusion regime this characteristic 
is the number of molecules in the power 2/3. 

One can see in fig. 8 three different curves and some solitary points. The 
solitary points correspond to the numerical simulation of the quasihomoge- 
neous case and three curves correspond to three models in the quasihomoge- 
neous case. 

The shortest spectrum is in the first model. The careful look shows that 
this line is doubled. This occurs because there the ideal variant of the first 
model is also drawn. This ideal variant corresponds to Wf ree /W to tai = 1 
in the subintegral function. The coincidence means that the main role in 
the first model is played by the relatively big droplets which were formed at 

W free = W tota i. 

The longest spectrum corresponds to the second model. This curve is 
very close to the intermediate curve which corresponds to the third model. 
The approximate coincidence of the second and the third models shows that 
both of them are valid and the role of the relatively big droplets here is the 
main one. 

One also see that the first model isn't too far from the real solution. This 
allows to present the rigorous estimates for the nucleation rate. Certainly 
the first model is the estimate of the real process from below. It gives the 
number of the droplets about 20 percents less than the numerical simulation. 

An estimate for the nucleation rate from above can be gotten by the 
following way. From the first model it follows that until t = 0.52 (practically 
this case is drawn in fig. 2) the rate of nucleation is near the ideal value 
and the deviation is less than 15 percents. So, one can say that the period 
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< t < 0.52 corresponds to the absence of overlapping (the first model is 
the estimate from above). Then one can consider the process where the total 
volume is exhausted only by the ER of the droplets appeared at < t < 0.52 
in a random manner. The distribution of the centers of ER of such droplets 
is evidently random. This model certainly gives the estimate from above for 
the nucleation process. The total number of the droplets is only 15 percents 
greater than the result of the numerical simulation. As a conclusion one can 
state that two suitable estimates from below and above are obtained. 

The proximity of the last estimate to the real solution justifies the sup- 
position that the main role in vapor consumption belongs to the droplets of 
the relatively big sizes appeared in the system practically free from the ER. 
This supposition can be also justified in the analytical manner. 

One can see that the second and the third models are rather close to the 
real solution but don't coincide with it. There are at least two reasons of 
such deviation. The first one is the presence of the strong correlations in 
the real system - if two ER overlap in some moment of time then the power 
of overlapping can only grow in time. It hasn't a random character as it is 
stated in the second and third models. 

This effect can be taken into account in a rather simple manner. It is 
sufficient to consider two spheres and calculate the power of overlapping as 
the function of the distance and time (it is the simple geometrical problem). 
Unfortunately the answer can be written only in a very complicate form. If 
we have two ER with radii R\ and R 2 correspondingly at a distance / between 
their centers and I > max{R 1 , R 2 ) then the volume of overlapping is 

2>kR\. n , . 2-kRI,^ n , , 

Vover = (1 - 2cOS(f 1 + COS </?i) H (1 - 2cOS(f 2 + COS lf 2 ) 

where 

-R\ + RI + 1 2 

COSifi = — " 



COS(f 2 



2R t l 
R\ + R\ + I 2 



2R 2 l 

Certainly, this result can not lead to the compact form of the balance equa- 
tion. It will be difficult to solve it analytically. 

The second reason is the " moderate overlapping" problem. This property 
means that actually there is an interaction through overlapping in ensemble of 
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several droplets. Earlier this property was extracted [O] in a case of a special 
effective length of ER. Now we see that this property has a rather general 
sense. The way to solve this problem proposed in [IT| is very complicate and 
leads to some uncertain relations. 

How one can overcome all these problems? Actually one has no need to 
do it. The simple numerical simulation takes into account all these effects 
and gives the universal solution. Really we need only one number - the total 
number of appeared droplets. It is known. Then one can forget about all 
mentioned difficulties. 

Now one can analyze the heterogeneous case explicitly. The suitable 
approximation is given by (|13|). One has to substitute instead of JV/, m the 
number of droplets given by the corresponding model. 

The relative error of approximation ([13]) is drawn in fig. 9 for the first 
model, in fig. 10 for the second model, in fig. 11 for the third model. It is 
rather small for all models. For the third model it is practically negligible. It 
is clear because the third model is based on the free volume approximation. 
Then the source of error can be found only in exhaustion of centers without 
any influence on the growing of ER. Without this influence the power of 
exhaustion can not essentially act on the process of nucleation. 

One can fulfill the same analysis for the numerical simulation. In fig. 12 
the relative error of ( |T3"D for numerical simulation is drawn. Here in (|T3|) 
the value Nh om = 504.8 from the numerical simulation is used. The result 
is compared with the computer simulation of heterogeneous condensation. 
This simulation is rather simple. One can take the procedure for the quasi- 
homogeneous case but put the center of the new droplet with probability 
rj/rjtot- Every time when this point will be out of ER we shall reduce rj as 
77 — >■ 77 — 1. 

One can see that the relative error is very small. We don't use the aver- 
aging over (this is the reason why there is no smooth curve) many attempts 
in order to see that the error of (|T^) has the scale of the mean square error 
of numerical simulation^. So, there is no need to use the more sophisticated 
approach. 

The solution of the problem now is completed. 

The generalization for the conditions of dynamic type is absolutely analo- 
gous to M. The convergence due to avalanche consumption is more weak and 
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For the system with 500 droplets. 
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one has to use instead of approximation ( |T3"D more sophisticated procedure 
described in |J. Certainly, the universal constants have to be calculated by 
the numerical simulation with profiles. 

The generalization for the arbitrary regime of the droplets growth can 
be done absolutely analogous to |§. This generalization is based on the 
similarity of the functional forms obtained here and in the AA. This similarity 



lies in the base of universality formulated in [ 1(| , JTTJ] . 

One can see that the theory of condensation with profiles taken into ac- 
count presents the picture which is absolutely different from the AA. Never- 
theless in many situations the result of experiment coincides with the result 
of the AA. One has to explain this coincidence. It is rather formal one. 

Really, in any experiment it is more convenient to have a small system 
and to get many droplets. The rate of nucleation has to be taken as a rather 
high one. So, the supersaturation is relatively high and parameter a -1 isn't 
a real small parameter of a theory^. Then as it is shown in the AA 
gives the write qualitative result despite the wrong base of consideration. 
The reason lies in the fact that at small a the main quantity of a substance 
is lying in the tail of profile. The tail of profile is rather thin and can be 



be taken into account by the AA. In ['1 1 the correction term for the AA at 
small a is also presented. 

The last important feature to mention is the movement of the embryos 
boundaries. This problem is widely discussed in the determination of the rate 
of regular growth for the supercritical embryos. In different systems the effect 
of influence of the boundary movement on the rate of growth is different . We 
have to note that in theory presented here the rate of the embryos growth is 
an external value which is assumed to be known^. 

Another problem is the adequate account of the effect of boundary move- 
ment in construction of the ER. One can see that in the first part of publica- 
tion it is already shown that the effect of the boundary movement is small. 
Here we shall present the abstract arguments for such conclusion. 

To use the thermodynamic approach the initial power of the mother phase 
metastability has to be relatively small. Together with the Maxwell's rule 
it leads to the following final result of the phase transition: " Only relatively 



24 It isn't necessary for consideration here, but it has to be small due to the possibility 
of thermodynamic description of the critical embryo. 
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It is really known practically for all systems. 



36 



small part of the system volume is occupied by the new phase." It isn't 
in contradiction with the property that all volume is occupied by ER. So, 
the final state of the system is the practically saturated mother phase in 
practically all volume of the system and the small volume (distributed over 
all system) occupied by a new phase. 

As the result one can see that the process of a substance consumption 
(extraction) leads to the saturation in a volume relatively large in comparison 
with a volume of the new phase embryo. 

The mother phase can not be undersaturated (then there will be the dis- 
appearance of embryos). So, as long as even the saturated phase has to 
be spread on large distances 1$ one can state that the growth of embryo 
produces the perturbation over the large distances I > Iq. To have an inter- 
ruption (relative interruption in comparison with an ideal nucleation rate) 
of a new phase formation one needs a very small reduction of the power of 
metastabilityQ. So, this reduction is attained at the distances V > I which 
are very large in comparison with the embryos linear size. Then one can use 
the point source approximation as it was done in the first part and forget 
about the boundaries movement. The negligible character of the boundary 
movement is proven now for all possible systems. 

The heat extraction and account of all other intensive parameters of de- 
scription can be done analogously to flT0|. 
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Fig. 1 
The form of f ((3). 
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Fig. 2 

The cross section of the system at t = 0.5 for the quasihomogeneous 

situation. 
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Fig. 3 

The cross section of the system att — 1 for the quasihomogeneous situation. 
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Fig. 4 

The cross section of the system at t = 1.5 for the quasihomogeneous system. 
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Fig. 5 

The cross section of the system at t = 1.5 for r\ tot = 50. 
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Fig. 6 

The cross section of the system at t = 3 for i] tot = 50. 
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Fig. 7 

The cross section of the system at t = 6 for i] tot = 50. 
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Comparison of different models in the quasihomogeneous situation. 
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Fig. 9 

Relative error of the quasihomogeneous approach in the first model. 



47 




Fig. 10 

Relative error of the quasihomogeneous approach in the second model. 
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Fig. 11 

Relative error of the quasihomogeneous approach in the third model. 
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Fig. 12 

Relative error of the quasihomogeneous approach in the universal 

simulation. 
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